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Abstract 

We present a systematic linear analysis of the perturbations in- 
duced by cosmic strings in cold dark matter. We calculate the power 
spectrum and find that the strings produce a great deal of power on 
small scales. We show that the perturbations on interesting scales 
are the result of many uncorrelated string motions, which indicates a 
much more “gaussian” distribution them was previously supposed. 

Cosmic strings are widely considered to be possible seeds for galaxies 
and larger scale structure (see, for example, these pioneering papers [1,2, 

3, 4, 5, 6, 7, 8, 9] and the numerous papers which have followed). The idea is 
that a network of strings formed as defects in a symmetry breaking phase 
transition in the early universe. Such a network is expected to rapidly 
approach a “scaling solution” in which most of the properties of the network 
are independent of the initial details. As the network evolves, the strings 
attract other matter gravitationally. In this way perturbations are seeded 
which can grow, via gravitational collapse, into galaxies and larger scale 

structure. _ t 

Veeraraghavan and Stebbins [10] have developed a formalism with which 
to study linear perturbations induced in the surrounding matter by sources 
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such as strings. We have applied this formalism to a numerical realiza- 
tion of a cosmic string network for the case of a flat universe with cold 
dark matter (CDM). We describe the results in terms of a simple model 
which can be easily extrapolated beyond the dynamical limitations of the 
numerical work, and adapted to different pictures of the string network. A 
complete discussion of this work will appear in [11]. This letter reports the 
main results, which differ in major ways from earlier results based on more 
heuristic approaches. 

The degree to which string networks have been shown to exhibit scaling 
behavior on cosmic timescales is still controversial [12]. For this work, 
however, we will assume that a realistic string network is well represented 
by the “scaling solution”. The scaling solution provides a very convenient 
way to think about a cosmic string network. At any given time the strings 
can be divided up into long strings and loops. Most of the properties of the 
long strings can be described in terms of a single scale £ defined by 

PL = f 2 , (1) 

where pi is the energy density in long string and p is the mass per unit 
length of the string. (Except where stated otherwise, we use geometrical 
units, where c = G — 1.) The long strings describe random walks of step 
size £, and also have a mean separation of 0(£), The strings move around 
at relativistic speeds, but the general properties of the network at two 
different times are related by simply rescaling £ appropriately. The degree 
to which the long string is straight on scales smaller than £ is a subject 
of current debate. One possibility is that small scale structure contributes 
to a renormalization of /z, but does not affect the gross scaling properties. 
In this case one can use the renormalized \l in Eq 1 so that £ reflects the 
large-scale spatial properties of the network. 

The time evolution of the string network is a process of equilibration, 
whereby the long strings are continually chopped up into the statistically 
favored loops (whose size is typically smaller than £). As a result pi de- 
creases, and £ grows. The expansion of the universe is described by the 
growth of the scale factor, a(t ). We consider the standard flat Friedmann- 
Robertson-Walker universe which has an intial radiation dominated equa- 
tion of state (^4> a(f) ot f ly/2 ) and which eventually becomes matter domi- 
nated (^i> a(t) oc t 2 / 3 ). We consider strings which are produced early in the 


radiation era. The canonical scaling picture dictates that, well into the ra- 
diation or matter eras, any initial conditions approach the scaling solution 
with ( <x R h ( Rh = a/a), but not necessarily with the same constant of 
proportionality. During the radiation-matter transition £/Rh interpolates 
smoothly between these two constant values. 

The quantity relevant to the perturbation theory is 6 ( 2 ) = (pd{x) ” 
Pd)/Pd, where pd{x) is the density of the dark matter today, and the overbar 
represents the spatial average. The (linear) perturbation induced in the 
dark matter may be written as the convolution of the string “source density 
with a suitable Green function, integrated over all time. The source density 
is written in terms of the string stress-energy density: ©+(*) = ©oo(z) + 
©jj(a;). In Fourier space 

St = Sl + 4tt(1 + Z eq ) f°° T(k; f?')©+(*b v'W ( 2 ) 

* k Jm 


where T is T 2 C from [ 10 ]. The effect of perturbations other than those due 
to string motions on the dark matter today is represented by Si. This 
term is responsible for “compensation”, in which string perturbations are 
canceled by other perturbations on sufficiently large scales, resulting in no 
net perturbation on scales too large to be affected by causal dynamics. 
The issue of compensation is a crucial one in these calculations and it is 
discussed in [ 11 ]. There we conclude that the important effects of Si, will 
be well approximated if the substitution T -* (1 + ^) _2 T is made and Si 
is dropped explicitly from Eq 2. The scale on which compensation sets in 

is given by k c . ... 

Thus modified, Eq 2 can be squared and averaged over directions in 

order to study P(k ), the power spectrum: 


(2n) 3 P(k)SW(k-k') = 16tt 2 (1 +z eq ) 2 f [ T(fc;7ji)T(fc';77 2 ) (©+(*, »7i)0 + (fc', 172 )) dr h d V 2 


The picture can be simplified by noting that although there are two inte- 
grals over conformal time, 77, the only significant contributions come from 
times when 771 and 772 are reasonably close. That is because the strings 
configurations are uncorrelated when separated by a sufficiently long time. 
For large enough values of t] x — r] 2 , (0 + (fc,T7i)©+(fc,q2)) will be negligibly 
small. As long as the correlation time of the strings is small compared with 
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the time over which T(k;r/) varies, we can model P(k) by 

roc 

P(k) = lftrV(l + **)’ / WijijWijijW (4) 

JT)i 

where 

(2tv) 3 P(k\ r)')6(k - k’) ^ J ^0 + (fc,T7')0 + (fc',7/")^dT7". (5) 

The “structure function”, T{k\7]), represents the power spectrum of the 
coherent string motions. In order to represent a scaling string network, 
we write F{k\r]) = F(£fc). (For some of our calculations T does vary by 
factors 0(1) over the string correlation time, but we do not expect this to 
introduce major errors, other than factors of 0(1) in the normalization.) 

One of our important results is that the contributions to the power spec- 
trum come from string motions at many different times (as one integrates 
Eq 4). A correct calculation of the full spectrum requires knowledge of the 
string network far beyond the dynamic range of any string simulation. This 
is where scaling is a big help. In order to calculate the power spectrum, all 
we really need from the simulations is the correct form for F(£k). This only 
requires an understanding of the network over one coherence time (but one 
still must resolve a range of spatial scales). 

We have chosen a general form for F( £k) which is motivated by cosmic 
string physics [13]: 


" j(lTiW)' (6 > 

We show in [11] that the model fits the Albrecht and Turok (AT) string 
network very welL Since the large scale fluctuations in the string density 
are just the result of uncorrelated fluctuations on the scale the structure 
function takes on a constant (or “white noise”) form for A >> £ (where 
A = 2ir/k), The coherent motions of the long strings on scales 0(£) and 
smaller trace out sheet-like wakes, and these result in a part of F(£k) 
which goes as fc“ 2 , which is characteristic of two dimensional surfaces. The 
structure function can be further modified to include the effects of loops, 
but we find in [11] that the resulting effects are negligible on the scales of 
interest. 
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For our calculations, the nature of the cosmic strings is completely spec- 
ified once one fixes the compensation scale k c and the function Civ)- We 
consider three different choices for these quantities. The AT parameters 
are k c = 4w /p and £( r ?) given by Eq 3.7 in [14]. The “ATO” (AT old) 
parameters are k c = 4x jr] and £(77) — Rh I n/IO- The ATO parameters rep- 
resent the old picture of cosmic strings (as extrapolated from [15] in the 
simplest way), and we use them to illustrate how little our results depend 
on the changing understanding of string networks. The BB-AS (Bennett 
and Bouchet [16] and Allen and Shellard [17]) parameters are k c = 2tt/t) 
and £(tj) = rj. The BB-AS simulations tend to show large values ([Rh We 
have intentionally exaggerated their results in this direction by choosing 
£ = rj (= H = “the horizon size”), in order to test the robustness of our 
conclusions. (For additional information on the string simulations see the 

contributions on that subject in [18].) 

We stress that £ represent the large scale geometrical features of the 
network so the renormalized fi appears in 6, incorporating the smaller scale 
features of the string. This renormalization must represent the contribu- 
tions of the small scale features to the effective linear source denstiy, as 
discussed in [19]. 

We plot (47r)(A/(27r))- 3 P(A) in Figure 1 (AT = short-dashed, BB-AS 
= solid, ATO = long-dashed), along with the same curve for perturbations 
from inflation in CDM [20] (dot-dashed). (We have used for the hubble 
constant today: H 0 = 1/R H = 100h km sec" 1 Afpc" 1 ) All the curves have 
the standard normalization (meant to reproduce the observed galaxy-galaxy 
correlation function) where the rms hmjm in an R — 8 Mpc sphere is unity 
(see Eq 8). This means that p is different for each curve (see the table 
below). The slope of all the curves is A" 4 for large A, because the scaling 
properties of the strings induce the same behavior in P( A) as is produced 
by the scaling properties of the inflationary perturbations. The transition 
between radiation and matter domination exposes the differences between 
strings and inflation. Smaller wavelengths, which enter the horizon in the 
radiation era have much more power in the string spectrum, compared with 
inflation spectrum. 

Stebbins et. al [21] have shown that wakes produced around the time of 
equal density in matter and radiation (t eq ) have the largest surface density. 
This fact has been widely interpreted to mean that the dominant long string 
contributions to the perturbations on many scales come at around t eq . We 
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have checked this directly by starting the time integration at t eq . The result 
is much less power on small scales. For example, the BB-AS power is about 
a factor of ten smaller at A — IMpc when the integration starts at t eq . The 
dominant contributions to small scales are produced before t eq . 

We fit the curves in Fig 1 to the form 


47r k 3 P(k) 


a X fil47rk 4 f 1 \ 2 

1 + (a 2 &) + (a 3 fc) 2 + (tt 4 ^) 3 \1 + (a 5 fc) 2 / 


( 7 ) 


The last factor (involving a 5 ) represents a deviation from A -4 behavior 
at very large wavelengths. This deviation occurs because the very largest 
scales have just entered the horizon and have yet to receive their full compli- 
ment of perturbations. In Table 1 we present values of the a’s which provide 
a fit to within 10% of the calculated curves: The values of /x r , which enters 



ai 

a 2 

a 3 

a 4 

a 5 

Mr 

AT 

45.0 

2.30 

1.80 

0.520 

.00127 

0.178 x 10- 6 

BB-AS 

21.8 

7.57 

5.89 

1.93 

.000357 

0.661 x 10~ 6 

ATO 

37.0 

8.95 

4.80 

1.25 

.00112 

0.360 x 10“ 6 


Table 1: Fit parameters (a’s) for Eq 7 when k is in units of Mpc J ,The 
values of p r are fixed by the normalization of P(k). 


as a parameter in Eq 6, are determined by standard 8 Mpc normalization 
condition. The subscript w r” emphasizes that this is the renormalized //, 
which incorporates possible effects of the small scale wiggles. 

We turn now to a more detailed analysis of our results. The integration 
over conformal time in Eq 4 covers the entire history of a given scale. At 
early enough times, any given scale is far outside the Hubble radius and is 
unperturbed. Around the time when R# (which always grows faster than 
a(f )) catches up with the comoving scale A, the string perturbations become 
uncompensated and P(A) starts to receive significant contributions. The 
power continues to receive contributions at all later times, but perturbations 
produced later have less time to grow, and this can diminish their relative 
contribution to P( A) 
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It will be useful to define the average of any quantity A , by inserting 
a factor X in the integrand of Eq 4 and then dividing the result by P{ A). 
We write the resulting average (X)\. In Fig 2 we plot the average redshift, 
( z/z eq )x 5 versus wavelength. One can see that different scales get their main 
contributions at different times. No one time (such as z eq ) is singled out. 

Although the perturbation theory is performed in Fourier space, the 
string network is best thought of in position space. The link can be made 
by the standard formula for mean squared mass perturbation, Sm 2 jm , 
inside a sphere of radius R: 


6m 2 
m 2 


J w 2 (x)P{k)A'Kk 2 dk 


( 8 ) 


where x = Rh~ l k , and w(x), the “window function” is given by: 

, , 3(sin(x) — *cos(®)) rQ >, 

w i x ) = ^ • [V) 

We thus define the position space average of any quantity A , by inserting 
a factor of X in the integrand of 8 and dividing by 6m 2 /m 2 . The result is 

labeled { X)r . . „ 

Figure 3 shows {N^r where N ( is the “number of £-times ”, 

B sfwj, - tjj ». (10) 

(?/A 

The string motions are coherent on a timescale At/ ~ so Fig 3 shows that 
most scales are perturbed significantly over more than one string coherence 
time (although for the BB-AS strings N ( is 0(1)). The relative magnitudes 
of the three different curves corresponds to the different relative sizes of 
The characteristic time over which a given scale is perturbed is similar 
in the three cases. (The slower coherent string velocities reported in [18] 
for the AS simulations may suggest a longer coherence time. However, for 
their large values of £, scaling actually requires that a segment of length £ 
loose a large fraction its length to loops in just one £-time . We feel that 
the presence of such a disruption on the timescale £ indicates that £ is a 
good measure of the coherence time even for the BB-AS strings.) 

If A > £, the scale A is being perturbed by many separate string motions 
occurring on the smaller scale £. A box of size A contains (A/ £) 3 boxes of size 
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£, so one can estimate the number of separate string motions (or “ coherent 
string wakes”) which contribute to a given scale as N c = ^{(A/^) 3 )*, and 
we plot (N c )r in Fig 4. One can see that many different wakes contribute 
to the perturbation on a given scale. 

If a string produces a perturbation at a given redshift, z , in the matter 
era, that perturbation has grown by a factor of z by today. Thus the 
earliest possible perturbations incur the largest growth. One reason we 
get large values of N c is that the perturbations for A > £ are produced 
earlier, and thus contribute more to the total perturbation. (Even though 

= 0(1) for the BB-AS strings, the contributions come mainly when 
A > £, resulting in large values of N c .) In contrast, perturbations produced 
in the radiation era do not grow much until the onset of matter domination, 
so the “growth factor” flattens out. The fact that strings significantly 
perturb the surrounding matter on scales smaller than the horizon combined 
with the flatter growth factor causes the dominant contributions to come for 
A < £ on scales which entered the horizon far enough back in the radiation 
era. This is the reason why N c gets smaller for small R. It is also the reason 
why the string spectrum has more power on small scales than the inflation 
spectrum* If the perturbations on a given comoving scale are dominated 
by the effects of the strings just when that scale crosses the horizon, then a 
spectrum similar to the inflation spectrum results [5]. The fact that in the 
radiation era the strings add significant contributions well after the horizon 
crossing boosts the power on small scales relative to inflation. 

As originally pointed out by Vachaspati [22], it is tempting to think that 
individual string wakes could be responsible for the sheet-like structures in 
the distribution of galaxies observed on scales of around 50 h~ l Mpc [23]. For 
this effect to arise, one needs a single coherent string motion to produce 
the dominant perturbation on a range of scales, so that the single wake 
dominates inside many spheres of different radii which contain the wake. 
Figure 2 is not encouraging in this regard, since it tells us that the average 
contribution to the perturbation on each scale comes at different times. 

If one centers a sphere of radius R on a planar string wake then, using 
[21], the perturbation today (in linear theory) is given by 6m jm — |A R~ l , 
where A is the distance out to which matter has been accreted. (A depends 
on the time when the wake was produced, and is proportional to p T s.) For 
a wake to produce a distinctive effect, Smfm around the wake must stand 
out above the rms Smfm on a range of scales. The straight lines in Figure 5 
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show the fjm/m around individual wakes. The curves show the rms 6m /m. 
For each set of string parameters, we calculate two string wake lines. One 
corresponds to when £ = bOMpc and the other is chosen so that A takes 
on its maximum value (A = 1.5 x 10 6 p, r h~ 2 Mpc), which occurs at about 
z = 2 z eq (see [21]) and corresponds to different £’s for different string pa- 
rameters. When R is about half the mean separation of wakes (« £/2), the 
value of 6m jm around the wake should be no different than the contribu- 
tion to the rms value of 6m jm produced by those wakes. We terminate the 
single- wake lines at R = £/2. Depending on what fraction of the rms 6m jm 
a particular wake contributes, the endpoint of the corresponding line may 
lie close to or below the rms 6m jm curve. (The three £ = 50 Mpc lines may 
be identified because they all terminate at R = 25Mpc.) We should note 
that since the single wake lines do not include the effects of compensation, 
they are actually upper bounds. 

Figure 5 shows that even though the perturbation around a given wake 
increases with decreasing R, the rms 6m jm is increasing at a similar rate. 
Since the slope of the rms 6m jm is close to that of the single wake 6m jm, 
the individual wake does not stand out except on scales orders of magnitude 
smaller than 50 h~ x Mpc. On the larger scales, it gets strong competition 
from the total 6m jm, which includes contributions from different wakes 
at different times. (It is possible that when loops are properly accounted 
for the single wake lines will not rise above the rms 6m jm on any scale, 
but if they do, that fact could be relevant to galaxy formation.) To the 
extent that many uncorrelated string motions contribute to the pertur- 
bation spectrum, the perturbations are more “gaussian” than previously 
suspected. We discuss this issue further in [11]. 

In conclusion, our picture of perturbations from cosmic strings in CDM 
has been significantly changed by a more systematic analysis. We have 
shown that the perturbation spectrum is the result of many uncorrelated 
string motions occurring at different times. The resulting spectrum has 
much more power on small scales, as compared with inflation, and also as 
compared with earlier pictures of perturbations from strings. The contri- 
butions of many string motions make it unlikely that distinctive stringy 
features will show up strongly in the matter distribution today. Our con- 
clusions assume a scaling string network, but otherwise depend very little 
on the string network parameters. 

Although the relationship between the linear perturbation spectrum 
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and the astrophysical observations is not well understood, we note that, 
based on current thinking, the string spectrum does not look good in this 
regard. The extra power in the string spectrum on small scales feeds into 
the normalization (that is, it contributes to the galaxy-galaxy correlation 
function). As a result the amplitude on large scales is reduced, as compared 
with inflation (see Fig 1). To the extent that the inflationary CDM spec- 
trum is in trouble for lack of power on large scales (as argued, for example, 
in [24]), one would expect the string induced CDM spectrum to have even 
greater problems. 

We thank Robert Brandenberger for useful conversations. This work 
was supported in part by the DOE and the NASA (grant NAGW-1340) at 
Fermilab. 
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Figure Captions 


Figure 1: The logarithm of 47r(A/27r)- 3 P(A) vs. log(A) for inflation (dot- 
dashed), AT strings (short-dashed), BB-AS strings (solid) and ATO strings 
(long-dashed). We use h = O 0 = 1 
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Figure 2: The the average redshift at which each scale is perturbed 

({ z/z eq )\ ) vs. A for AT strings (short-dashed), BB-AS strings (solid) and 
ATO strings (long-dotted). 


Figure 3: The number of string coherence times ((N{)r) vs. R for 

AT strings (short-dashed), BB-AS strings (solid) and ATO strings (long- 
dashed). 


Figure 4: The number of contributing coherent string motions ((N c )r) vs. 
R for AT strings (short-dashed), BB-AS strings (solid) and ATO strings 
(long-dashed). 


Figure 5: The value of 6m jm in a ball of radius iZ, vs. R for AT strings 
(short-dashed), BB-AS strings (solid) and ATO strings (long-dashed). The 
curves give the rms £m/ra, and the straight lines give 6m jm in a ball 
centered on a single wake. We use h = = 1 
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